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ABSTRACT 

We find self-similar solutions to describe the interaction of spherically sym- 
metric ejecta expanding at relativistic speeds with an ambient medium having 
a power law density distribution. Using this solution, the time evolution of the 
Lorentz factor of the outer shock is derived as a function of the explosion energy, 
the mass of the ejecta, and parameters for the ambient medium. These solutions 
are an ultra-relativistic version of the solutions for the circumstellar interaction 
of supernova ejecta obtained by Chevalier and extensions of the relativistic blast 
wave solutions of Blandford Sz Mckee. 

Subject headings: gamma rays: bursts — hydrodynamics — relativity — shock 
waves — supernovae: general 



1. Introduction 

Self-similar solutions for fluid flows involving shock waves propagating at relativistic 
speeds have been presented in the literature (e.g. Blandford & McKee 1976; Hidalgo & 
Mendoza 2005; Nakayama & Shigeyama 2005). Most of them are ultra-relativistic versions of 
the already known non-relativistic self-similar solutions. Blandford & McKee (1976) derived 
solutions for an ultra-relativistic spherical blast wave enclosed by a strong shock, which is 
initiated by a release of a large amount of energy within a small volume. If the released 
energy is not so large, the resultant non-relativistic flow is described by the well known 
solution discovered independently by Sedov (1959) and Taylor (1950). The implosion of a 
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spherical strong shock wave had been solved independently by Landau and Stanyukovich 
(see e.g. Zel'dovich & Raizer 2002) and by Guderly (1942) in non-relativistic cases. Hidalgo 
& Mendoza (2005) derived self-similar solutions for the same problem in the relativistic limit. 
Nakayama & Shigeyama (2005) discovered solutions for the emergence of an ultra-relativistic 
shock wave at the surface of a star. When the shock wave is non-relativistic the same problem 
had been posed by Gandel'man & Frank-Kamenetskii (1956) and solved by Sakurai (1960). 

Here we present a self-similar solution for the interaction of spherical ejecta expanding 
at ultra-relativistic speeds with the ambient medium. When the ejecta expand at non- 
relativistic speeds, self-similar solutions are known to exist (Chevalier 1982). He discovered 
solutions for the flows resulting from a collision of a freely expanding spherical matter with 
the ambient medium. In the relativistic limit, self-similar solutions describing the shocked 
ambient media with power law density distributions were already obtained by Blandford & 
McKee (1976) with the shock Lorentz factor evolving with time t as t~ m l 2 . The constant 
m can be determined when the nature of explosion is specified. This paper discusses self- 
similar solutions to describe shocked ejecta as a result of a collision with the ambient medium. 
This solution will be relevant to extremely high energy supernovae originated from compact 
progenitors with small masses. 

The density distribution of the ejecta used here is suggested by analytical solutions 
(Nakayama & Shigeyama 2005) for the breakout of an ultra-relativistic shock wave from a 
hydrostatic atmosphere and its non-relativistic versions (Sakurai 1960; Matzner & McKee 
1999). The outer layers of the resultant flow in a non-relativistic solution expand at ve- 
locities greater than the speed of light. They have a density distribution that agrees with 
numerical calculations for the same problem using a special relativistic hydrodynamical code 
(Nakamura & Shigeyama 2004). Thus we use the results as the inner boundary conditions 
of the present self-similar solutions. We will describe the relation of our solution with results 
obtained by Blandford & McKee (1976) in some detail. 

In the next section, the density distribution of freely expanding ejecta as a result of 
shock breakout is discussed. Self-similar flows in the shocked ejecta and ambient medium 
are obtained in §3. The behavior of these two flows at the contact surface is discussed in §4. 
Using the obtained solution, the shock propagation in the ambient medium is discussed in 
§5. The results are summarized in §6. 
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2. Relativistic ejecta from shock breakout 

The velocity v c of freely expanding ejecta at a distance r from the explosion site is 
written as v c = r/tai a time t after the explosion. The corresponding Lorentz factor is 

c 



7 e = \j yj\ — v\. We set the speed of light equal to unity throughout this paper. In the 
relativistic limit, self-similar solutions for the emergence of a shock wave at the stellar surface 
(Nakayama & Shigeyama 2005) suggest that the density in the outer layer of ejecta has a 
power law distribution in terms of the Lorentz factor as, 

where p e is the density of ejecta measured in the co-moving frame and p' e in the fixed frame, 
that is, p' c = p c 7 (hereafter densities denoted with primes mean that they are measured in 
the fixed frame). Here the non-dimensional parameter n can be estimated to be ~ 1.1 if the 
stellar envelope is in radiative equilibrium and the adiabatic index is 4/3. The constant a has 
the dimension of a mass. On the other hand, a self-similar solution for the non-relativistic 
shock breakout in a plane-parallel geometry (Sakurai 1960; Fields et al. 2002) results in 
the density distribution in the same form with n ~ 2.6. Though this self-similar solution 
assumes non-relativistic flow, the resultant density distribution agrees with that obtained 
by numerical calculations solving special relativistic hydrodynamics equations (Nakamura 
& Shigeyama 2004) even for large Lorentz factors 7 C > 10. The transition between these 
solutions occurs when the pressure to the mass density P/p at the shock front becomes of 
the order of unity. Thus the outermost layers should be governed by the solution in the 
relativistic limit. 



3. Flow in the shocked region 

Flows located between two shock fronts are considered by connecting two self-similar 
solutions. Suppose that the shock front in the ambient medium propagates at a speed V\ 
and the shock front in the ejecta at V 2 . If we define a = 7/(7 — 1), where 7 is the adiabatic 
index, the flow in each of the shocked regions is described by the following equations. 

, „ « d In '"y dP 9 dP , , 

dln(P a -V) a dr 2 v 
dt r 2 dr ' 
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where the pressure is denoted by P, the density p, the velocity v, the corresponding Lorentz 
factor 7, and the density in the fixed frame p' = p'-f. In the shocked ambient medium, 
the ultra-relativistic equation of state is assumed. Thus 7 = 4/3 and the energy density 
e per unit mass is written as e = 3P/ p by neglecting the rest mass energy density. These 
approximations lead to the same equations as presented in Blandford & McKee (1976). 
On the other hand, it is found that the rest mass energy should not be ignored in the 
shocked ejecta to have a compression shock. Otherwise the jump condition would result in 
an expansion shock for m < 1. Then the energy density becomes e = Pj {(7 — l)p} + 1 in 
this region and the term including the density in the left hand side of equation (2) is not 
neglected. 



3.1. Jump conditions at the shock front in the ejecta 

The density, pressure, and velocity have discontinuities at the shock front. Relations to 
connect these quantities at both sides of the shock front are given here. The mass flux pw M 
and the energy momentum tensor T^ u change across the shock propagating at the velocity 
of V 2 according to the following formulae, 

[pu^ = 0, (5) 
[T^K = 0, (6) 

where [F] = F c — F 2 . The subscript e refers to values in the freely expanding ejecta at radius 
R 2 defined as 

R 2 (t) = [ drV 2 (r), (7) 
Jo 

while the subscript 2 refers to values in the shocked matter at the same radius R 2 . The 
variables n M , u^, and T^ v have been defined as 

TV ee T 2 (-V 2 , 1, 0, 0), u» = 7(1, v, 0, 0), (8) 
T» v = phu^u u + Prf, (9) 

with the metric rf v = diag(— 1, 1, 1, 1). Here T 2 denotes the Lorentz factor of the shock 
and is assumed to evolve as a function of time t as 

= Ar m (10) 

with constants A and m. The enthalpy h is denoted by other thermodynamic variables as 
h = e + Pj p. Equation (5) yields 



Pc7c(^ c - V 2 ) = p 2 -f 2 (v 2 -V 2 ), 



(11) 
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Equations (6) can be rewritten as 

p el 2 c (v c - V 2 ) = p 2 h 2l \{v 2 - V 2 ) + P 2 V 2 , (12) 

p c 1 2 c v c (v c - V 2 ) = p 2 h 2 -f 2 v 2 (v 2 - V 2 ) + P 2 . (13) 

Equations (11)-(13) together with the relation 7 e = y/m + 1T 2 at r = R 2 lead to the following 
equation 

ax 3 + (a- 2)\/m+ lx 2 -(a- 2)x - ay/m + l = 0, (14) 

where x = r y 2 /T 2 . This equation has the only one positive solution between 1 and \/m + 1 
for m > — 1. If this solution is denoted as x = y/g, the Lorentz factor, density, and pressure 
at the front in the shocked matter will be derived as 

72 = ?ri, (is) 

P2 = 7 ^T7 TTPe7e, 16 

(m + l)(q - 1) 



These relations are used to obtain the boundary conditions for flows in the shocked ejecta. If 
the constant m is smaller than unity as is usually the case, the compression is so weak that 
the thermal energy cannot dominate the rest mass energy in the shocked ejecta. Therefore we 
cannot use the ultra-relativistic equation of state there while the ultra-relativistic equation 
of state is always a good approximation in the shocked ambient medium as long as the shock 
Lorentz factor is much larger than unity. 



3.2. Self-similar flow in the shocked ejecta 

It is convenient to define self-similar variables for the pressure F(£), the Lorentz factor 
G(£), and the density H(£) in the shocked ejecta as 



p = tt^-„ M-r]FK>. US) 



a(q- 1) + 2 V V m + !, 

7 2 = (19) 

The similarity variable ^ is defined as 

e = {l + 2(m+l)l1}(l-^), (21) 
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following Blandford & McKee (1976). 

To obtain equations governing the self-similar evolution, the transformation of coordi- 
nates from (r, t) to (T 2 , £) is performed using the following relations of derivatives; 

<|: = -{i + 2 ('" + i) r I}^. (23) 

d d , _ x / 1 A d , niS 

-™7^7TT2+(™+ 1 ) 7^7- (24) 



dint dlnT 2 v ' \qG S J <9£ 

Then we obtain from equations (2)-(4), 

2( m+ i)(i + gGO^^ - (m + l)(a + 7)(1 - qG^)^^- = {mn - (a + J)m - 6}gG, (25) 

<i In F d In G 

2(a-l)(m+l)(l-gG£)— jt a(m+l)(l+gGf)— nr- = {(l-a)mn+am+2a-6}gG', (26) 

2(m + 1)(1 - gGO^TT^ - 2 ( m + !)^^ = ~( mn ~ m - 2 )?G, (27) 
where / = /(£) = p/P. The boundary conditions are given at the shock front as 

G(l) = F(l) = H(l) = 1. (28) 

The solution is obtained by numerically integrating equations (25)-(27) from the shock front 
£ = 1 to the contact discontinuity £ = £ c = l/qG(£ c ). 



3.3. Jump conditions at the shock front in the ambient medium 

The density, Lorentz factor, and pressure change across the shock wave propagating at 
the speed V\ = y/l — 1/T 2 according to the following relations, 

A = e£,* = £, Pl = Ml, (29) 

where the subscript 1 refers to values in the shocked fluid at the shock front, while the 
subscript i refers to values in the unshocked ambient medium at the shock front. The 
shock Lorentz factor is denoted by T. This T evolves with time in the same way as 1^, i.e., 
T 2 = Bt~ m . The density in the ambient medium is assumed to have a power law distribution 



-7- 



as p = br k where b is a constant. Thus self-similar variables for the pressure /(x), Lorentz 
factor g(x)i an d density h(x) are defined as 

P= 2 -^ff(x), (30) 

7 2 = y<?(x), (31) 

p' = 2 A r 2 / i ( X ). (32) 
Here the similarity variable x nas been introduced as 

X={l + 2(m + l)r 2 }(l-^). (33) 



3.4. Self-similar flow in the shocked ambient medium 

The same procedure in §3.2 leads to equations governing a self-similar flow in the shocked 
ambient medium as 

d\nf A{2(m-l) + k}-(m + k-A)g X _ ^ 



gdx (m + 1) (4 - Sgx + g 2 X 2 ) 

ding (7m + 3A; — 4) — (m + 2)gx 
gdx (m + 1) (4 - 8gx + g 2 X 2 ) 
dlnh _ 2(9m + 5k-8)-2(5m + Ak-Q)gx+ (m + k - 2)g 2 x 2 
~gd~x~ ~ (m + l)(A-8gx + g 2 X 2 )(Z-gx) ' 

in the relativistic limit. The boundary conditions are 



(35) 
(36) 



/(l)=0(l) = fc(l) = l. (37) 
These are exactly the same equations as presented in Blandford & McKee (1976). 

4. Flow at the contact surface 

The density distribution has a discontinuity at the contact surface defined as £ = £ c and 
X — Xc where G(£ c )£c = 1/q and g(Xc)Xc — 2, while the velocity and pressure are continuous. 
The same velocity of the two flows at the contact surface indicates that the ratio of Lorentz 
factors of the two shock fronts should become 

T 2 B G(L) , x 
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according to equations (19) and (31). The continuous pressure at the contact surface requires 
that the pressures at both sides of the contact surface evolve with time in the same manner. 
Thus we obtain from equations (18) and (30) 



This indicates that k < 4 + n/2 because m > — 1. In addition, the pressure equilibrium 
yields a relation of constants that have been introduced in the preceding sections as 



This relation will be used to discuss the propagation of the outer shock in the ambient 
medium. 

A blast wave in a uniform ambient medium with power supply was treated by Blandford 
& McKee (1976). They assumed that the power is supplied by a stationary source located at 
r = and varies with time as a power law, L = Lot s . If s = (2n — 8)/(n + 2), then the flow 
in the shocked ambient medium evolves exactly in the same way as described in Blandford 
& McKee (1976). They also argued that the adiabatic impulsive solution is appropriate for 
the flow with m in the range of m > 3 or m < —1. Since m > — 1 in our problem, the 
condition m > 3 restricts the value of n to be less than —2k/3 through equation (39). Thus 
the values of n = 1.1 and 2.6 indicate a continuous supply of energy. 

When (n, k) = (1.1, 2), equation (39) yields m = 2/(n + 2) ~ 0.65 and then, from 
equation (14), q ~ 1.06. Performing numerical integration of equations (25)-(27) and (34)- 
(36) with these parameters, boundary conditions (Eqs. (28) and (37)), and 7 = 5/3 in the 
shocked ejecta, we find £ c ~ 0.95 and (/ c , g c , F c , G c ) ~ (0.18, 0.47, 2.1, 0.98). Here the 
subscript c of each variable indicates that the value is taken at the contact surface. Results 
are shown in Figure la. When (n, k) = (2.6, 2), the same procedure leads torn = 2/(n + 2) ~ 
0.43, q ~ 1.05, £ c ~ 0.97, and (/ c , g c , F c , G c ) ~ (0.32, 0.70, 2.0, 0.98). Results are shown in 
Figure lb. 



M l+n/2 



a 




(40) 



5. Propagation of shock in stellar wind 



Nakayama & Shigeyama (2005) investigated the evolution of an ultra-relativistic shock 
wave in a plane-parallel atmosphere and derived self-similar solutions. The resultant energy 
spectrum of the ejected matter with an explosion energy E ex and an ejected mass M e j for a 
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radiative stellar envelope can be expressed as 



5.75 



where M(> e) denotes the ejected mass that have particle energy per nucleon greater than 
e = m p (7 e — 1). Here the proportionality constant is not specified by Nakayama & Shigeyama 
(2005) because their solutions cannot trace the flow to the free expansion phase in which the 
equation of state is no longer ultra-relativistic. 

Using the relation dM = Anr 2 p c "i c dr we can convert the energy distribution given by 
equation (41) into the density distribution of the ejecta as a function of 7 C and t. 

Pc = a 7e -r 3 oc 5 ' 75 M e]1 ;^h-\ (42) 

This formula indicates a oc E^ 5 M~^ 4 ' 75 in equation (1), n = 1.10, and m = 0.645 from 
equation (39). Eliminating A in equations (38) and (40) leads to 



r 2 



3m(2g)"/ 2 ( & \ n/2 ( ^ _ \~^~ \ F c G n J 2 a 



2{a(q -1) + 2}V^ + V V V m + lj f c g^ 2 b 



2/(n+2) 

r m . (43) 



Substituting the expression for a into equation (43) with the exponents m and n obtained 
above, then we have 

T oc {a/b) 1/{n+2) t' m l 2 oc £^ 5 M7 L53 r a322 . (44) 

We can discuss the time evolution of the Lorentz factor of a shock in a steady stellar 
wind more definitely if the proportionality constant in equation (41) is specified. This is 
the case for Type Ic supernovae (SNe Ic) discussed by Fields et al. (2002) and Nakamura & 
Shigeyama (2004). According to the results of special- relativistic hydrodynamical calcula- 
tions by Nakamura & Shigeyama (2004) (see also Fields et al. 2002), the energy distribution 
of SN Ic ejecta is given by 



M(> g) = c ( E cx \ 3 - 4 / " 3 - 4 / e x -3-6 

M ej V 1051er g s / V 1 M &) VlOMeV/ ' 1 1 

where C is a non-dimensional constant (~ 2 x 10~ 4 ) slightly depending on stellar models. 
Thus the density distribution of the ejecta becomes 
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Accordingly the constant a is expressed in terms of the parameters of explosion as 

a = 2.26 x 10- 8 CM cj ( f cx V 'V -^O " . (47) 
J Vl0 51 ergs/ V 1 M 0/ 

On the other hand, the constant 6 is expressed in terms of the mass loss rate M and the 
wind velocity f w ind as 6 = M/47rt> wind . Substituting the expression for a obtained above and 
values for non-dimensional functions at the contact surface into equation (43), T is expressed 
in terms of physical quantities characterizing the explosion and ambient medium as 

r = Uxl0 (iw) (lo^J (nt) (t^J ■ (48) 

6. Summary 

We have found self-similar solutions for the collision of spherical ejecta freely expanding 
at relativistic speeds with the ambient medium when the ejecta have a density distribution 
in the form of equation (1) and the ambient medium has a power law density distribution 
in terms of the distance from the explosion site with the exponent —k. Solutions can be 
obtained as long as k < 4 + n/2. It is found that the rest mass energy is always comparable 
to the thermal energy in the shocked ejecta while the rest mass energy becomes negligible 
in the shocked ambient medium. 

Profiles of the flows in the shocked regions are presented in two cases where the ejecta 
are results of the shock breakout from a stellar surface embedded in the steady wind. The 
time evolution of the Lorentz factor of the outer shock front is derived as a function of 
physical parameters characterizing the explosion and the environment. 

We thank an anonymous referee for important remarks. This work has been partially 
supported by the grant in aid (16540213) of the Ministry of Education, Science, Culture, 
and Sports in Japan. 
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Fig. I. — (a): Profiles of the pressure P (solid line), the density p (dashed), and the gas 
Lorentz factor 7 (dash-dotted) in the interaction region between freely expanding ejecta with 
n = 1.10 and the ambient medium with a power low distribution (p oc r~ 2 ). The physical 
variables are normalized to their values at r = Ri where the subscript 1 refers to the value at 
the outer shock front. The radius at the contact surface is denoted by -R c . (b): Profiles of the 
same quantities as (a) but for freely expanding ejecta (E cx = 10 53 ergs and M C j = 13 M & ) 
at time t = 100 sec. 



